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ABSTRACT 


A  model  for  prediction  of  particle  radius  and  oxide  cap  size/shape  versus  time  for 
an  aluminum  particle  tracking  a  stream-tube  through  a  solid  rocket  motor  port 
has  been  developed.  Following  preliminary  calculations  leading  to  a  postulated 
flame  structure,  a  quasi-steady  model  to  predict  instantaneous  consumption  of 
aluminum  and  generation  of  condensed  oxide  (both  as  a  cap  on  the  aluminum 
particle  and  as  smoke)  as  a  function  of  instantaneous  particle  size,  ambient 
conditions,  and  cumulative  amount  of  oxide  in  the  cap  was  developed.  Finally, 
this  model  was  imbedded  into  a  framework  tracking  evolution  of  particle  size, 
oxide  cap  size,  and  ambient  conditions,  which  change  as  the  particle  travels 
along  a  stream-tube  consuming  oxidizer  and  releasing  heat.  Qualitative 
agreement  was  found  between  predictions  and  limited  observations. 


ABSTRACT 


A  model  for  prediction  of  particle  radius  and  oxide  cap  size/shape  versus  time  for 
an  aluminum  particle  tracking  a  stream-tube  through  a  solid  rocket  motor  port 
has  been  developed.  Following  preliminary  calculations  leading  to  a  postulated 
flame  structure,  a  quasi-steady  model  to  predict  instantaneous  consumption  of 
aluminum  and  generation  of  condensed  oxide  (both  as  a  cap  on  the  aluminum 
particle  and  as  smoke)  as  a  function  of  instantaneous  particle  size,  ambient 
conditions,  and  cumulative  amount  of  oxide  in  the  cap  was  developed.  Finally, 
this  model  was  imbedded  into  a  framework  tracking  evolution  of  particle  size, 
oxide  cap  size,  and  ambient  conditions,  which  change  as  the  particle  travels 
along  a  stream-tube  consuming  oxidizer  and  releasing  heat.  Qualitative 
agreement  was  found  between  predictions  and  limited  observations. 


BACKGROUND/PRELIMINARY  CALCULATIONS 

In  1995,  Brooks  and  Beckstead  [1]  published  a  paper  in  which  the  state-of-the-art 
regarding  aluminum  particle  combustion  modeling  and  data  was  well  reviewed. 
Subsequent  papers  by  Babuk  [2],  DesJardin  [3],  Fabignon  [4],  and  Bazyn  [5] 
added  further  information  in  this  area.  Due  to  word  limitations,  these  papers  will 
not  be  discussed  in  detail.  The  model  development  presented  in  the  next  section 
drew  on  all  of  these  papers  as  well  as  on  the  following  preliminary  calculations. 

(1)  Aluminum  partial  pressure  adjacent  to  the  particle  surface  was  determined  as 
a  function  of  surface  temperature  from  KP  =  PAi.g  (atm)  =  exp(-AG°/RTs)  where 
AG°  is  the  change  in  free  energy  accompanying  vaporization  as  tabulated  in  the 
JANAF  Thermochemical  Tables  [6], 

(2)  This  relationship  was  combined  with  a  simplified  B-Number  analysis  (as 
described  in  Turns’  textbook  [7])  to  estimate  the  surface  temperature:  this 
calculation  indicates  that  the  surface  temperature  is  far  below  the  boiling  point  of 
aluminum  at  normal  motor  operating  pressures,  counter  to  assumptions  made  in 
some  previous  models. 

(3)  From  several  sets  of  thermochemistry  calculations,  it  was  concluded  that  the 
predominant  AlxOy  specie  would  be  AI2O3  liquid  at  any  place  in  the  flame 
structure  where  the  oxygen/aluminum  atomic  ratio  is  high  enough  to  take  it  there. 

(4)  Further  free  energy  calculations  indicated  that  Babuk’s  [2]  claim  that  4AI(I)  + 
AI2O3  (I)  will  convert  to  3  AI2O  (g)  at  the  oxide-aluminum  interface  in  the  oxide 
cap  region  is  unlikely. 
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(5)  Additional  thermochemistry  calculations  indicated  that  likely  reactions  at  the 
aluminum  surface  are  Al  vaporization  and  4AIO  ->  Al203  (I)  +  Al20  (g). 

(6)  Still  more  thermochemistry  calculations  indicate  that  if  the  principal  oxidizers 
are  C02  and  H20  (typical  in  an  aluminized  solid  rocket  motor)  the  likely  reaction 
sequence  is  Al(g)  +  C02  ^  AI0/1/2AI202  +  CO  and  Al(g)  +  H20  ^  AI0/1/2AI202 
+H2  at  a  primary  flame  sheet  followed  by  further  oxidation  to  Al203  (I)  +  CO  +  H2 
at  an  outer  condensation  sheet. 

(7)  Finally,  thermochemistry  calculations  were  used  to  estimate  the  fraction  of 
AlxOx  which  will  appear  as  monomer  rather  than  dimer  as  a  function  of  flame 
sheet  temperature:  these  calculations  yielded  Mass  Fraction  Monomer  =  0.122  + 
0.000577(T-3500) 


MODEL  DEVELOPMENT 

A  sketch  of  the  postulated  quasi-steady  combustion  model  is  presented  as 
Figure  1 .  (As  will  be  discussed  later,  this  model  is  imbedded  in  another  code  in 
which  the  progressions  of  boundary  conditions  as  the  particle  travels  along  a 
given  stream-tube  are  updated.)  As  shown,  the  structure  is  composed  of  three 
zones  (Inner,  Middle,  and  Outer)  bounded  by  the  droplet  surface,  an 
infinitesimally  thin  flame  sheet  where  primary  oxidation  occurs,  and  another 
infinitesimally  thin  zone  where  final  oxidation  to  Al203  and  oxide  condensation 
occur.  It  is  postulated  that,  in  addition  to  vaporization  of  aluminum  from  the 
particle  surface,  AIO  diffusing  back  to  the  surface  from  the  primary  flame  sheet 
undergoes  collision-limited  reaction  with  the  surface  to  form  Al203  liquid,  which  is 
added  to  the  cap,  and  Al20(g)  (via  4AIO  ->  Al203(l)  +  Al20(g)),  with  the  aluminum 
vapor  and  Al20  then  traveling  outward  toward  the  primary  flame  sheet.  At  this 
flame  sheet,  it  is  assumed  that  the  aluminum  vapor  reacts  instantaneously  with 
C02  and  H20  (coming  from  the  outer  regions)  to  form  a  mixture  of  AIO  and  Al202, 
with  the  ratio  of  monomer  to  dimer  being  determined  by  the  flame  zone 
temperature,  as  discussed  earlier.  Based  on  the  preliminary  calculations,  it  is 
assumed  that  the  Al20  does  not  react  here  but  instead  travels  outward  along  with 
the  AI0/AI202  mixture  to  the  condensation  sheet  where  the  suboxides  react  with 
additional  carbon  dioxide  and  water  vapor  to  produce  Al203(l),  hydrogen,  and 
carbon  monoxide.  It  is  further  assumed  (flame  sheet  approximation)  that  there 
are  no  gas  phase  reactions  occurring  in  the  three  zones  permitting  standard 
development  of  algebraic  equations  describing  the  transport  of  species  and 
enthalpy  across  the  three  zones.  Fick’s  first  law,  with  “effective”  binary  diffusion 
coefficients  calculated  at  the  midpoint  temperature  of  each  zone  and  with  Lewis 
Number  =  1  is  used  in  development  of  these  equations. 


Total  enthalpy  (chemical  plus  sensible)  is  assumed  to  be  constant  across  each 
zone  (no  reactions)  with  the  values  of  the  total  enthalpy  in  each  zone  being 
determined  by  energy  balances  at  the  surface,  at  the  primary  flame  sheet,  and  at 
the  condensation  sheet,  accounting  for  source/sink  terms  associated  with 
reactions  at  those  locations.  In  addition,  the  outer  boundary  conditions,  which 
would  normally  be  specified  at  an  infinite  radius,  are  instead  applied  at  a  radius 
specified  by  the  spacing  of  the  particles  relative  to  one  another  in  the  stream- 
tube  flow.  An  “inert”  specie,  representing  all  the  free-stream  species  other  than 
C02,  H20,  CO,  and  CO  is  included  as  part  of  the  overall  gas  composition. 
Obviously,  for  a  quasi-steady  analysis  (no  accumulation  terms)  its  flux  through 
each  zone  will  be  zero,  as  will  fluxes  of  CO  and  H20  in  the  inner  zone,  where 
diffusion  and  convection  of  these  species  are  required  to  cancel  out.  The 
aluminum  oxide  smoke  produced  at  the  condensation  zone  as  a  ‘sink”  term  is 
assumed  to  disappear  from  the  system  though  its  thermal  effects  must  be 
included.  For  solution  of  this  problem,  the  condensation  temperature  has  to  be  a 
model  input,  based  on  the  previously  discussed  thermochemistry  considerations. 

Boundary  conditions  for  solution  of  these  equations  describing  the  aluminum 
burning  rate,  deposition  of  oxide  at  the  surface  (as  characterized  by  oxide  cap 
size)  and  production  of  oxide  smoke  at  the  condensation  “sink”  region  are 
brought  into  the  single  particle  analysis  from  the  wrap-around  code  in  which  this 
analysis  is  imbedded.  These  conditions  include  free-stream  temperature  (TAmb), 
free-stream  mass  fractions  of  C02,  H20,  CO,  CO,  and  “inert”  (Yj.amb), 
instantaneous  particle  radius  (rsurf),  fraction  of  surface  coated  by  oxide  (b), 
pressure  (P),  amount  of  aluminum  still  unburned  (Mai),  and  mass  of  oxide 
accumulated  in  the  cap  (Mox). 

In  line  with  the  assumption  of  infinitesimally  thin  reaction/condensation  sheets,  it 
is  assumed  that  the  mass  fractions  of  Al(g),  H20,  and  C02  at  the  flame  sheet  are 
zero  as  are  the  mass  fractions  of  Al20  and  AlxOx  at  the  condensation  sheet.  In 
addition,  the  ratio  of  AIO  to  Al202  is  assumed  to  be  frozen  at  the  flame  sheet 
radius  value  elsewhere.  The  remaining  unknowns  after  application  of  these 
assumptions/approximations  are  mass  fractions  of  species  and  mass  fluxes 
through  the  three  zones 

At  r  =  rsurf  :  Y Ai,s>  YaixOx,  s,  YCo,s-  YH2.s,  Y|NERT  S 

At  r=rFz  :  Yaio.fz,  YA1202,  fz,  Yai2o,fz,  Yco.fz,  Yh2,fz,  Yinert.fz 
At  r  =  rcz  :  Yh2.cz,  Yco.cz,  Yc02.cz,  Yh2o.cz,  Yinert.cz 

Inner  Zone  (rs  <  r  <  rFz) : 

Al , Inner  ’  AlxOx, Inner  ’  Al  20 , Inner  ’  Net , Inner 


Mid  Zone  (rFZ  <  r  <  rCz) 

^ Al2OMid’  ^ AlxOx  Mid’’  ,^II2OMid’’  ^ CO,Mid ’  ^CCX  Mid ’  ^NeiMid 

Outer  Zone  rCz  <  r  <  rAMB)  : 

^H2,OuteP  ^H20,0uteP  ^ CQOuteP  ^COJJuteP  ^NetjOutei 


Additional  unknowns  are  the  surface  and  flame-sheet  temperatures  (TsUrf  ,  TFz), 
total  enthalpy  flux  across  each  zone  (Hinner,  HMid,  and  Horn),  flame  sheet  radius 
(rFZ),  Condensation  Sheet  Radius  (rcz),  mass  rate  of  aluminum  consumption  at 
the  surface  (dMAL/dt),  mass  rate  of  oxide  deposition  at  the  surface  (dM0x/dt)  and 
the  mass  rate  of  production  of  oxide  smoke  at  the  condensation  radius  (itisink), 
for  a  total  of  42  unknowns. 


Accordingly,  42  equations  are  required  for  closure  of  the  problem,  given  the 
boundary  conditions  brought  in  at  each  instant  of  time  from  the  wrap-around 
stream-tracking  code.  Four  equations  are  provided  by  integration  of  the 
convective/diffusive  species  flux  equations  (in  the  absence  of  reaction)  across 
the  outer  zone: 


i¥ 


net ,  outer 


A7TpD{\-  p) 


wY  Y  —  rfif 

ft  I  -j  j  ^ _  1  net ,  outer  j  ,AMB  j, Outer 

Kllrcz~llrAU y  — ft 

'*net , outer  J  j ,CZ  j , outer 


(1-4) 


j  =  CO,  co2,  h2,  h2o 


Similarly,  six  additional  equations  are  provided  by  integration  of  the 
convective/diffusive  transport  equations  across  the  middle  zone: 


net  Mid 


4  npD{\  -  P) 


/I  /„  _l/„  \  _  1n  "*-net,Mid1  j,CZ  j  Mid 

V '  'FZ  l'rcz)-m  ft  y  _ft 

net  Mid 1  j,FZ  "  *  j  Mid 


(5-10) 


j  =  CO,  co2,  h2,  h2o,  aixox,  ai2o 


Five  more  equations  are  produced  through  integration  of  the  convective/diffusive 
transport  equations  across  the  inner  zone: 


net, Inner 
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v  -  M 
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net, Inner  j,Surf  j, Inner 


j  =  CO,  h2,  aixox  ,  ai2o,  ai 


where: 


for  each  zone 


(16-18) 


In  addition,  Y|nert  at  the  surface,  flame  sheet,  and  condensation  sheet  are 
obtained  from: 


2>.  =  i.o 


for  each  location 


(19-21) 


Integration  of  differential  equations  for  total  (chemical  plus  sensible)  enthalpy 
across  each  of  the  three  zones  (Inner,  Middle,  and  Outer  yields  three  equations 
relating  temperatures  at  r  =  rSurf,  r  =  rFz,  and  r  =  rCz  to  the  enthalpy  fluxes  in  each 
zone  (Hp, Inner,  Hf,mid,  and  HF, Outer)  Of  the  form: 
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iKj,Outefp,j 


47rpD{\-/3)cp 


(\/rcz-\/rAmh)  =  \n 


H F, Outer 

(Z flute fi p ,  j ) TCZ  +  Z ^jaj  ~  H F, Outer 


(22) 
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(23) 
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j, Inner  pj 
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(24) 


where  the  aj’s  are  (h°  —cpT°)j  values  obtained  from  curve  fits  for  each  species 

and  the  HF’s  are  obtained  from  enthalpy  balances  at  the  surface,  flame  sheet, 
and  condensation  sheet  as: 


H  =0593¥  h  +M  h 

1 JF ,Inner  7  jaA  lfix  Jnner  A  /203  (l),TSutf  T  niA!JiineflAKl)TSulf 

^  h\  Mid  ^ F, Inner 

H  =  H  -  rtf  h 

11  F, Outer  11  F, Mid  Sink  n Al 2O30\TCZ 

Next,  there  are  several  stoichiometric  relationships  to  be  invoked: 
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Finally: 
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(41) 


(from  collision  rate  theory,  with  v  being  the  input  fraction  of  collisions  of  AIO  at 
the  surface  resulting  in  reaction),  and 


(42) 


P(Yaim/27)/(£  Yjs /MW  j)  =  exp(  12 .43  -  34680  /Ts„rf) 

OUTLINE  OF  SOLUTION  PROCEDURE 

First,  values  for  the  ambient  temperature  (values  at  r  =  tamb)  and  the  ambient 
mass  fractions  of  CO,  C02,  H2,  H20,  and  inert  (all  other  species  being  lumped  as 
inert,  with  straightforward  calculation  of  an  average  molecular  weight  for  that 
inert)  were  imported  along  with  current  particle  radius  (rSurf),  pressure  (P),  and 
fraction  of  the  particle  surface  covered  with  oxide  (b)  from  the  stream-tube 
tracking  wrap-around  program  described  in  the  next  section.  Given  these  values, 
the  approach  taken  to  solving  this  set  of  equations  consisted  of  a  double-nested 
trial-and-error  loop  in  which  dMA|/dt  and  dM0x/dt  were  first  guessed  and  used  to 
sequentially  calculate  other  parameters  until  Yco2  at  the  flame  sheet  was 
calculated.  If  this  answer  was  not  0.0,  the  guess  on  dlVWdt  was  adjusted  and 
this  inner  loop  was  repeated  until  the  C02  mass  fraction  at  the  surface  did  equal 
zero.  As  part  of  this  loop  a  value  for  the  condensation  radius  was  calculated. 
Next,  the  remaining  equations  were  used  to  arrive  at  a  second  value  for  the 
condensation  radius  based  on  the  enthalpy  balances.  If  the  two  values  of 
condensation  radius  did  not  agree,  a  new  value  of  dMox/dt  was  chosen,  the 
entire  inner  loop  was  reconverged,  and  another  comparison  of  the  condensation 
radius  calculated  in  the  inner  loop  with  a  new  value  calculated  from  the  enthalpy 
balances  was  made.  This  overall  procedure  was  continued  until  both  loops  were 
converged.  During  all  of  these  loops,  updates  on  surface  temperature,  flame 
sheet  temperature,  and  midpoint  temperature  values  for  each  zone  were 
calculated  and  used  to  update  values  on  the  density-diffusivity  products  in  each 
zone  and  the  mean  heat  capacities  in  each  zone.  In  addition,  the  average 
molecular  weights  of  the  gas  mixtures  in  each  zone  were  continuously  updated. 
More  details  on  this  approach  can  be  obtained  from  the  author.  Hard-wired 
values  of  the  various  cpj’s,  a/s,  and  values  of  species  enthalpies  at  a  reference 
temperature  were  built  into  the  computer  code. 

Only  on  the  first  call  from  the  wrap-around  stream-tube  tracking  code  is  this 
process  somewhat  time-consuming  if  bad  initial  guesses  are  made.  On 
subsequent  calls,  made  at  small  time  intervals,  the  changes  in  ambient 
conditions,  particle  radius,  and  fraction  coating  of  oxide  are  quite  small  and  the 
final  converged  solution  for  the  previous  time  step  is  most  adequate  for  rapid 
convergence  in  the  present  time  step. 

WRAP-AROUND  STREAM-TRACKING  CODE 

The  model  described  above  is  imbedded  in  a  stream-tracking  code  which  keeps 
track  of  all  of  ambient  and  particle  parameters  as  the  particle  travels  through  the 
rocket  motor,  heating  (releasing  heat)  and  eating  (consuming  oxidizer)  as  it  goes. 
Inputs  to  this  code  include  the  initial  particle  size,  initial  oxide  thickness, 
pressure,  stagnant  film  radius  (based  on  particle  spacing),  and  several  initial 
guesses  for  parameters  being  calculated  in  the  imbedded  code  at  zero  time. 


Included  in  the  code  are  curve  fits  for  ambient  species  mass  fractions  and 
ambient  temperature  as  functions  of  the  cumulative  percentage  of  aluminum 
burned  for  the  propellant  formulation  of  interest.  Curve  fits  for  condensed 
aluminum  and  aluminum  oxide  densities  as  functions  of  temperature  are  also 
included.  In  addition,  a  rather  simple  geometry-based  algorithm  for  the 
calculation  of  the  fractional  surface  coating  by  the  molten  oxide  cap  as  a  function 
of  the  instantaneous  masses  of  aluminum  and  oxide  cap  is  included. 

The  densities  of  liquid  aluminum  and  aluminum  oxide  are  built  in  as: 

pAi  (gm/cm3)  =  2.364[1 .0  -  0.00002268(T  -  933)]  (43) 

PA1203  (gm/cm3  =  2.998[1 .0  -  0.00003731  (T  -  2327)]  (44) 

where  T  is  in  degrees-Kelvin,  and  are  used  to  relate  the  aluminum  particle  radius 
to  the  aluminum  mass  and  the  oxide  cap  volume  to  the  cap  oxide  mass. 

The  geometry  of  the  oxide  cap,  which  leads  to  a  relationship  between  the 
fractional  aluminum  surface  area  coverage  and  the  volume  (and  thus  mass)  of 
oxide  in  the  cap  and  volume  (mass)  of  the  aluminum  sphere  was  calculated  by 
assuming  that  the  intersection  of  the  Al  and  Al203  is  represented  by  the 
intersection  of  two  spheres  at  a  90-degree  angle  (probably  somewhat  of  an 
oversimplification)  as  depicted  in  Figure  2.  For  a  90°  intersection,  it  can  be 
shown  that  the  distance  between  the  centers  of  the  spheres  (d)  is  related  to  the 
aluminum  droplet  radius  (R)  and  the  oxide  lobe  radius  (r)  by: 

d  =  ^R2  +r2  (45) 

and  that  the  intersection  of  the  spheres  is  a  circle  with  radius: 

a  =  rR/^R2  +  r2  (46) 

It  can  also  be  shown  that  the  volume  common  to  both  spheres  is: 

V  =  n(R  +  r -  df  (d1  +  2dr-  3r2  +  2dR+ 6  rR  -2>R2)H2d  (47 ) 

and  that  the  volume  of  the  oxide  cap  is: 

VOLOX=(4/3)xr}-V  (48) 

Since  the  cap  radius  (r)  is  not  initially  known,  but  VOLOX  is  known  via  the 
cumulative  oxide  cap  mass  accumulation,  an  iterative  procedure  is  used  to  find  a 


cap  radius  for  a  given  value  of  the  aluminum  sphere  radius  which  gives  the 
correct  value  for  VOLOX.  The  linear  angle  a  (see  Figure  2)  is  then  calculated  as: 

a  =  arcsin(  a  / R)  (49) 

and  the  solid  angle  as: 

SA  -  2k{\  -  cos  a)  (50) 

The  fraction  of  the  aluminum  sphere  radius  covered  with  oxide  is  finally 
calculated  as: 

/?  =  SA  /4  k  =  (1  -  cos  a)  12  (51) 

As  mentioned  earlier,  as  the  aluminum  particle  burns,  it  releases  heat  and 
consumes  oxidizer  from  its  surroundings.  To  treat  these  changes  in  the  ambient 
conditions,  correlations  of  the  various  species  mass  fractions,  temperature,  and 
molecular  weight  of  the  “inert”  were  generated  as  functions  of  the  cumulative 
percentage  of  aluminum  burned,  assuming  no  mixing  with  adjacent  stream  tubes 
or  diffusion/conduction  of  species/heat  along  the  stream  tube.  For  the  particular 
formulation  of  interest  (70/18/12  AP/AI/HTPB),  the  following  correlations  (with  X 
being  the  cumulative  weight  percent  aluminum  burned)  were  built  into  the  code: 


Y„  =  0.008 +  0.00022X 

(52) 

YHp  =  0.270- 0.001 3 4X 

(53) 

Yco=0.\  80+  0.00 1 60X 

(54) 

YCOi  =  0. 1 74  -  0.00243  IX  +  0.0000 1 04 7X2 

(55) 

Yinert  =  0.3  68 + 0.00 1 95 IX-  0.0000 1 047T2 

(56) 

MW INERT  =  33.28  -0.0185X-0 .000242 X2 

(57) 

TAmb  =  2550.  + 1 0.687X-  0.03 1 88X2  for  x  <  70 

T.imb  =  31 10.+ 14.666(^-70.) 


for  X  >  70 


(58) 


The  overall  operation  of  this  code  is  outlined  below: 


(1 )  Calculate  initial  mass  of  aluminum  and  oxide  cap.  Set  Time  =  0 

(2)  Calculate  initial  oxide  lobe  radius  and  fraction  of  aluminum  surface 
covered  ((3) 

(3)  Calculate  initial  ambient  temperature  and  composition. 

(4)  Call  aluminum  combustion  subroutine,  getting  back  mass  consumption 
rate  of  aluminum  and  generation  rates  of  oxide  cap  and  oxide  smoke. 

(5)  Take  a  small  time  step  and  calculate  masses  of  aluminum  burned  and 
oxides  generated  over  that  step. 

(6)  Calculate  new  value  of  aluminum  mass  and,  from  that  and  the  initial  mass, 
cumulative  percentage  burned. 

(7)  Calculate  new  values  of  aluminum  and  oxide  cap  volumes 

(8)  Using  Equations  52-58,  calculate  new  ambient  condition  values. 

(9)  Using  the  geometric  analysis  outlined  earlier,  calculate  a  new  value  for 
fraction  of  the  aluminum  particle  surface  area  blocked  by  the  oxide  cap. 

(10)  Call  the  aluminum  combustion  subroutine  again  and  continue  with  Steps 
5-10  ad  infinitum.  Track  cumulative  amount  of  aluminum  remaining, 
percent  burned,  cumulative  oxide  smoke  production,  and  cumulative 
amount  of  oxide  in  the  cap. 


SAMPLE  PROBLEM  DEFINITION  AND  OUTPUT 

The  sample  problem  entails  a  70/18/12  AP/Aluminum/HTPB  formulation  burning 
at  68  atmospheres,  producing  100-micron  radius  aluminum  particles  with  a  one- 
micron  thick  initial  oxide  coating,  which  is  assumed  to  immediately  draw  up  into  a 
lobe  on  one  side  of  the  particle.  The  input  aluminum  oxide  condensation 
temperature  was  4000K  and  the  input  reaction  efficiency  of  AI0/AI202  gas 
molecules  striking  the  virgin  aluminum  surface  was  unity.  The  input  effective 
ambient  radius  around  each  particle  was  set  based  on  the  aluminum  particle 
spacing  in  the  flow  through  the  rocket  motor  (for  this  case,  approximately  10 
times  the  initial  particle  radius.)  Results  are  presented  in  Figures  3-8. 

The  percentage  (mass)  of  the  particle  consumed  is  plotted  against  time  from 
ignition  out  to  approximately  120  milliseconds  in  Figure  3.  The  percent 
consumed  increases  monotonically  with  time,  with  a  constantly  decreasing  slope 


of  the  curve.  Thus,  approximately  55%  of  the  mass  has  been  consumed  at  40 
msec,  70%  at  60  msec,  90%  at  100  msec,  and  98%  at  120  msec  (where  the 
calculation  was  terminated). 

In  Figure  4,  the  aluminum  particle  radius  and  the  radius  of  the  oxide  lobe  are  also 
plotted  against  time.  The  aluminum  radius  is  seen  to  decrease  (slightly  less  than 
linearly)  with  time,  while  the  oxide  cap  radius  increases  monotonically  with  time, 
rapidly  at  first  but  at  a  decreasing  rate  with  increasing  time.  Beyond 
approximately  60  mesc,  the  oxide  cap  lobe  radius  is  larger  than  the  aluminum 
particle  radius. 

Since  in  classical  diffusion-limited  droplet  combustion,  a  d2-law  applies,  with 
droplet  radius  squared  decreasing  linearly  with  time,  the  square  of  the  predicted 
particle  radius  was  plotted  against  time  in  Figure  5.  A  significant  deviation  from 
d-squared  behavior  is  observed  with  the  absolute  value  of  the  slope  of  the  curve 
decreasing  significantly  with  time.  This  behavior  is  largely  a  result  of  the  fact  that 
the  fraction  of  the  aluminum  surface  blocked  by  the  oxide  cap  increases  with  time 
as  shown  in  Figure  6. 

As  discussed  earlier,  part  of  the  aluminum  oxide  appears  in  the  form  of  smoke  at 
the  condensation  zone,  while  part  deposits  on  an  accumulating  cap  at  the 
aluminum  particle  surface.  The  split  is,  of  course  of  interest  due  to  the  very 
different  hydrodynamic  behaviors  of  the  two  fractions.  Accordingly  fraction  of  the 
oxide  appearing  in  the  cap  is  plotted  against  time  in  Fig.  7 — this  is  seen  to 
asymptote  to  about  20  percent.  Finally,  in  Figure  7,  sketches  of  the  geometries 
of  the  aluminum  and  the  oxide  cap  are  presented  on  a  common  scale:  as  may  be 
seen,  the  final  oxide  cap  radius  is  predicted  to  be  approximately  70%  of  the 
original  aluminum  particle  radius. 
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PRELIMINARY  CALCULATIONS 


Aluminum  vapor  pressure  versus  temperature  (chart  on  Slide  3) 

Estimation  of  burning  aluminum  surface  temperature  and  mass 
fraction  of  Al  adjacent  to  the  surface  using  a  very  simplified  B- 
Number  analysis.  Indicates  surface  temperature  considerably 
below  boiling  point  and  surface  mass  fraction  Al  <0.10. 

Calculation  of  pseudo-boiling  point  of  Al203.  Ran  several  different 
types  of  thermochemistry  calculations,  concluding  that  the 
predominant  AlxOy  species  will  be  Al203  liquid  any  place  in  the 
flame  zone  where  the  oxygen/aluminum  atomic  ratio  is  high 
enough  to  take  it  there. 

Examination  of  Babuk’s  claim  that  4AI(I)  +  Al203  (I)  will  convert  to 
3AI20(g)  at  the  interface  between  Al  and  oxide  in  the  cap  region. 
Free  energy  calculations  indicate  that  this  will  not  occur  at 
reasonable  values  of  the  surface  temperature. 
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PRELIMINARY 


Definiti^AlHgiPrU)i#^Vflhl:^(N3ction(s)  at  the  aluminum  particle 
surface.  A  series  of  thermochemistry  calculations  at  specified 
temperatures  for  various  ratios  of  aluminum,  oxygen,  and  argon  mixtures 
indicate  that  at  the  surface  these  reactions  will  be  aluminum  vaporization 
and  4AIO— >  Al203(l)  +  Al20(g) 


Additional  free  energy  calculations  indicated  that  for  an  environment 
where  the  principal  oxidizers  are  C02  and  H20,  the  likely  sequence  of 
reactions  involves  Al(g)  +  C02  -->  AI0/1/2AI202  +CO  and  Al(g)  +  H20 
>  AI0/1/2AI202  +  H2  at  a  primary  flame  sheet  followed  by  further 
oxidation  of  Al20  +  AIO/1/2  Al202  with  C02  or  H20  to  produce  Al203(l) 
+  CO  and  H2  at  an  outer  condensation  sheet. 


Examination  of  the  fraction  of  AlxOx  which  will  appear  as  monomer 
versus  dimer  as  a  function  of  flame  sheet  temperature  (little  effect  on 
species  transport  but  a  significant  effect  on  flame  temperature).  Results 
of  these  calculations  are  shown  in  Slide  4  and  built  into  the  aluminum 
combustion  code  described  later. 
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FIGURE  1.  ALUMINUM  VAPOR  PRESSURE  CURVE 
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MASS  FRACTION  AL202  DISSOCIATED  TO  ALO 


TEMPERATURE  (I>EG-K) 


FIGURE  2.  A10-A1202  SPLIT  VS  FLAME 
ZONE  TEMPERATURE 
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MAJOR  ASSUMPTIONS  AND  APPROXIMATIONS 


■Infinitesimally  thin  reaction  sheets  (infinitely  fast  reactions)  at  primary  flame 
sheet  and  condensation  sheet.  AIO  reaction  rate  at  surface  equal  to  collision 
rate  multiplied  by  input  collision  efficiency 

■  Quasi-steady-state 

■Fick’s  Law  applies,  with  all  Lewis  Numbers  =  1 

■Constant  p,  D,  cp,  MWm,x,  etc.  across  each  zone,  evaluated  at  midpoint  T 

■No  gas-phase  reactions  in  zones  between  the  sheets:  permits  analytical 
integration  of  flux  equations  to  algebraic  forms 

■Al20  does  not  oxidize  further  at  flame  sheet,  but  reacts  at  condensation  sheet 
■Magical  “sink”  region  at  condensation  zone  for  aluminum  oxide  smoke 
■AIO  monomer/dimer  split  determined  by  temperature  at  the  flame  sheet 
■No  radiative  heat  transfer 

■  Input  value  for  T  at  the  condensation  sheet--fortunately  weak  dependence 


Distribution  A:  Public  Release;  distribution  unlimited 


Liquid  Aluminum  Droplet 


Surface 


Primary  Flame 
Sheet 


Inner  Zone  (No  Rx) 


Al(g) 


Al20(g)- 


■AICKg) 


Stagnant  Inert,  CO,  H2 


Net  Flux 


r  =  rsuif 

Al  valorizes  and 
4AI  ->  At203fl)  +  Al20(g> 


Condensation 

Sheet 


H2  1 

CO- 


CO2 

H2O 


Stagnant  Inert 


Net  Flux 


r  =rFZ 

Al  +  C02  ->0*10^  +  CO 
AI  +  t^O  ->(AIOh  +  h2 


Outer  Boundary 


Middle  Zone  (No  Rx) 

AI2O - 

AIO/AI2O2 - 


Outer  Zone  (No  Rx) 

-  CO2 

-« -  h2o 


CO 


H2O 


Stagnant  Inert 


Net  Flux 


'  =rcz 

Al20  +  2002  ->AI2Q3(]>  +  2C0 
Al20  +  2H20  ->AI2Q3a)  +  2^2 
2AI0+  002 ->AI2Q3a)  +  CO 
2AI0  +  H2O  >AI2Q3d)+  H2 
AI2P2  +CO2  >AI2Q3d)  +  CO 
Al202  +  H20  ->AI2Q3a)  +  H2 


'  =  'AMB 


Figure  3.  Sketch  of  Aluminum  Combustion  Model 
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GENERAL  OUTLINE  OF  SOLUTION  PROCEDURE 


■End  up  with  42  equations  in  42  unknowns:  rrij’s,  mNet’s,  and  enthalpy 
fluxes  across  each  zone,  Y-s  at  surface,  primary  flame  sheet,  and 
condensation  sheet,  surface  T,  primary  flame  sheet  temperature,  radial 
locations  of  the  primary  flame  sheet  and  the  condensation  sheet,  rate 
of  changes  in  aluminum  mass  and  oxide  cap  mass,  and  rate  of 
production  of  oxide  smoke 

■For  given  instantaneous  ambient  conditions,  particle  size,  and  oxide 
cap  size,  guess  the  rate  of  consumption  of  the  aluminum  droplet  and 
accumulation  of  the  oxide  cap,  and  use  a  double-nested  trial-and-error 
approach  to  home  in  the  correct  values 

■Update  property  values  as  part  of  each  loop 
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WRAP-AROUND  STREAM-TUBE  TRACKING  CODE 


■  Next,  the  code  described  on  the  previous  charts  was  imbedded  as  a  called  subroutine  in  a 
“wrap-around  code”  which  keeps  track  of  the  “heating  and  eating”  of  the  environment 
(increased  ambient  temperature  due  to  heat  released  by  the  Al  combustion  and  decreased 
oxidizer  level  due  to  C02  and  H20  consumption  in  the  combustion). 

■  First,  external  calculations  of  temperature  and  composition  of  the  stream-tube  as  a  function  of 
the  percent  aluminum  burned  for  the  given  propellant  formulation  and  pressure  are  made, 
and  curve  fits  of  temperature,  mass  fraction  C02,  mass  fraction  H20,  etc  as  functions  of  the 
percent  aluminum  burned  are  developed  and  built  into  the  “wrap-around  code”  which  keeps 
track  of  progress  of  all  variables  along  the  stream-tube.  In  addition,  curve-fits  of  aluminum 
and  aluminum  oxide  densities  vs  temperature  are  built  into  this  code. 

■Also  included  in  the  “wrap-around  code”  is  an  algorithm  for  calculation  of  the  changing 
geometry  of  the  aluminum  and  its  oxide  cap  as  mass  is  removed  from  the  aluminum  and 
oxide  cap  mass  increases  with  travel  along  the  stream-tube.  This  algorithm  permits 
updating  of  the  fractional  aluminum  surface  area  coated  by  the  cap  (P)  as  Al  combustion 
progresses.  This  algorithm  involves  assumption  that  the  Al  droplet  is  spherical  and  that  the 
oxide  lobe  is  a  sector  of  a  sphere  which  intersects  the  Al  surface  at  a  90  degree  angle  as 
depicted  in  the  next  chart  (possibly  an  oversimplification  which  will  have  to  later  be  modified). 
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Figure  4.  Sketch  of  Postulated  AI-AI2O3  Lobe  Geometry 
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SAMPLE  PROBLEM  DEFINITION  AND  OUTPUT 


■The  sample  problem  analyzed  entails  a  70/18/12  AP/AI/HTPB 
propellant  burning  at  68  atm  (1000  psia) 

■The  chosen  aluminum  particle  radius  was  100  microns 

■An  initial  uniform  oxide  coating  thickness  of  1  micron  was  selected, 
with  assumption  that  it  is  drawn  up  into  a  cap  at  the  time  of  ignition 

■The  input  oxide  smoke  condensation  temperature  was  chosen  to  be 
4000  degrees-K 

■It  was  also  assumed  that  every  collision  of  AIO  with  the  aluminum 
surface  results  in  reaction 

■The  input  effective  film  radius  around  each  particle  was  based  on  the 
aluminum  particle  spacing  in  the  flow  through  the  rocket  motor, 
approximately  10  times  the  initial  particle  radius  for  this  case 

■Results  are  presented  on  the  remaining  charts 
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Figure  5.  Percentage  Aluminum  Burned  Versus  Time 
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Figure  7.  Radius  of  Al  Sphere  and  Oxide  Cap  Lobe  Vs  Time 
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Figure  8.  Square  of  Radius  of  Al  Sphere  Vs  Time 
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Figure  9.  Fraction  of  Al  Surface  Area  Blocked  By  Oxide  Cap  Vs  Time 


Distribution  A:  Public  Release;  distribution  unlimited 


Fraction  of  Oxide  Product  in  Cap 


Time  (msec) 

Figure  10.  Percentage  of  Oxide  Product  Appearing  in  Cap  Vs  Time 
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FIGURE  1 1 .  Evolution  of  AI-AI2O3  Shapes 
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